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The possibility of making precise predictions for the Casimir force is essential for the theoretical 
interpretation of current precision experiments on the thermal Casimir effect with metallic plates, 
especially for sub-micron separations. For this purpose it is necessary to estimate very accurately 
the dielectric function of a conductor along the imaginary frequency axis. This task is complicated 
in the case of ohmic conductors, because optical data do not usually extend to sufficiently low fre- 
quencies to permit an accurate evaluation of the standard Kramers-Kronig integral used to compute 
e(i^). By making important improvements in the results of a previous paper by the author, it is 
shown that this difhculty can be resolved by considering suitable weighted dispersions relations, 
which strongly suppress the contribution of low frequencies. The weighted dispersion formulae pre- 
sented in this paper permit to estimate accurately the dielectric function of ohmic conductors for 
imaginary frequencies, on the basis of optical data extending from the IR to the UV, with no need 
of uncontrolled data extrapolations towards zero frequency that are instead necessary with stan- 
dard Kramers-Kronig relations. Applications to several sets of data for gold films are presented to 
demonstrate viability of the new dispersion formulae. 

PACS numbers: 05.30.-d, 77.22. Ch, 12.20.Ds 
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I. INTRODUCTION 



The Casimir effect, and more in general dispersion 
forces, are the current object of much theoretical and 
experimental interest. The reasons for the continuing in- 
terest in this well-established field are numerous, and are 
connected with both fundamental and applied science. 
For an updated overview of the subject, the reader is 
referred to recent reviews [H-Hl- 

After fifty years of slow progress, the field of Casimir 
physics received a strong push in the last decade, as a re- 
sult of a new wave of experiments 4:| which succeeded for 
the first time in measuring the tiny Casimir force with un- 
precedented precision. The new measurements provided 
a definitive confirmation of this phenomenon, but at the 
same time experiments utilizing metallic surfaces (which 
are used in the majority of current experiments) have 
raised unexpected theoretical puzzles, for the so-called 
thermal Casimir force. The problem is whether in the 
computation of the Casimir force at finite temperature 
conduction electrons in the plates should be described 
as a collisionless plasma, similarly to what is done in 
infrared optics, or if relaxation processes should be in- 
cluded The two approaches have been dubbed 
in the Casimir literature as plasma and Drude prescrip- 
tions, respectively. The problem is rather subtle indeed, 
with non trivial implications in thermodynamics and sta- 
tistical physics that are not yet fully understood. In fact. 



* Bimonte@na. infn. i1%l 



it has been found that inclusion of relaxation effects in 
the idealized case of two plates with no imperfections 
leads to violation of the Nernst theorem [8], while no 
such violation occurs when any amount of impurities is 
present [9]. It has been shown however that neglect of 
relaxation processes leads to contradiction with the Bohr- 
van Leeuwen theorem of classical statistical physics [lo| . 
From an experimental point of view, the resolution of 
the problem requires observation of the thermal Casimir 
force between two metallic plates, because the two alter- 
native theoretical pictures give different predictions for 
its magnitude. The experimental situation is perplexing, 
and at the moment of this writing there appears to be 
a striking contradiction between the results of different 
experiments by different groups. A series of experiments 
using microtorsional oscillators by the Purdue group [TJ] , 
providing at this time the most precise measurements of 
the Casimir force between metallic bodies in the sepa- 
ration range from 160 nm to 750 nm, were shown to be 
in agreement with the plasma prescription and to rule 
out the Drude prescription. Also a large distance torsion 
balance experiment |12| in which the Casimir force was 
measured in the separation range from 0.48 to 6.5 /im, 
obtained results that are in agreement with the assump- 
tion of ideal metal plates, and are in contradiction with 
the Drude model. On the contrary, a new experiment 
[l3| in which the Casimir force between a spherical lens 
and a flat plate was measured in the range from 700 nm 
to 7 fiui by a torsional balance, was claimed to be fully 
consistent with the Drude prescription, and to rule out 
the plasma one. This state of things calls for a careful 
investigation, to see if the contradiction can be resolved. 
It is well known by now f3| that the theoretical in- 
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terpretation of precision Casimir experiments is a difS- 
cult problem, because several possible sources of error 
must be considered, that include for example consider- 
ation of unavoidable imperfections of the surfaces, like 
roughness and other uncontrolled geometric distortions 
, the possible presence of residual electrostatic forces 
etc. We shall not be concerned with these issues here, 
and for a thorough discussion of the problem of theory- 
experiment comparison in Casimir ex peri ments, we ad- 
dress the reader to the recent work [TBlJlGj . In this paper 
we shall focuss our attention on the influence of the opti- 
cal properties of the the plates, described by their dielec- 
tric function e(u;), which by far constitutes the dominant 
factor affecting the magnitude of the Casimir force, espe- 
cially at sub-micron separations. The effect of the optical 
properties of the plates on the Casimir force must be esti- 
mated very precisely, if one seeks to discriminate between 
alternative prescriptions for the thermal Casimir effect, 
because the predicted magnitudes of the force resulting 
from different approaches differ by just a few per cent at 
separations below one micron, which is the range inves- 
tigated in the experiments by the Purdue group. 

According to the theory of dispersion forces [l| in or- 
der to compute the Casimir force between two dielec- 
tric (non-magnetic) plates of any shape at temperature 
T one needs to know the electric permittivities e(i^n) of 
the materials constituting the plates, evaluated at cer- 
tain temperature-dependent discrete imaginary frequen- 
cies = 2TrnkBT/h, where n is a non- negative integer, 
known as Matsubara frequencies. Getting accurate val- 
ues for the dielectric function e(i^) is not easy, in par- 
ticular when ohmic conductors are considered. This is 
so because the imaginary-frequency dielectric function 
e(i^) cannot be directly measured by any experiment, and 
therefore it has to be computed using Kramers-Kronig 
(KK) dispersion relations, starting from optical data of 
the involved materials. We shall specifically consider the 
case of gold, which is the material of interest in the re- 
cent experiments probing the thermal effect. Until a few 
years ago (T7l - [l9j it was customary to compute e(i£ ) by 
using the data for gold tabulated in the handbook [20]. 
However, in the most recent experiments the new prac- 
tice has been adopted of computing the quantities e{i£,n) 
on the basis of optical data measured directly on the 
samples used in the Casimir measurements. The neces- 
sity of doing this was first stressed in [2lj and it stems 
from the fact that optical properties of thin polycrys- 
talline gold films, typically used in Casimir experiments, 
are much dependent on the used deposition method. It 
has been claimed 0, Hsj that possible sample-to-sample 
variations of the optical properties may engender signifi- 
cant variations in the Casimir force, possibly larger than 
10 per cent at submicron separations. While the large 
differences in the Casimir force reported in (22| probably 
do not arise for carefully prepared gold films like those 
used in the experiments it seems correct to state 
that the optical properties of the used samples should 
be carefully characterized, if per cent precision in the 



theory-experiment comparison is aimed at. 

A big practical problem remains however because, due 
to the characteristic 1 /uj singularity displayed by the di- 
electric function of ohmic conductors, the standard KK 
integral used to compute e(i^) (Eq. ^ below) receives 
a very large contribution from frequencies below the IR, 
for which it is very difficult to collect optical data for 
thin metallic films. The remedy to this problem, ini- 
tially adopted in and followed afterwards by all au- 
thors, consists in extrapolating the available data below 
the minimum accessible frequency Wmin, in order to esti- 
mate the KK integral in the interval < lu < Wmin where 
optical data are not directly available. The extrapolation 
is usually done by means of the simple Drudc formula for 
the dielectric function of ohmic conductors. It is impor- 
tant to observe that this contribution is usually very large 
(it easily accounts for more than fifty per cent of the KK 
integral) and therefore the values of the Casimir force 
that are obtained by this procedure are very sensitive to 
experimental uncertainties in the values of the Drude pa- 
rameters. For example, in Ref. [22f| it was estimated that 
experimental uncertainties in the Drude parameters may 
imply by themselves an uncertainty in the Casimir force 
as large as 1 % at 100 nm. 

In our opinion, the reliability of the whole procedure 
described above can be questioned, because in reality 
there is no way to know how accurately the Drude model 
describes the dielectric function of gold films, in the wide 
range of frequencies for which it is used in Casimir com- 
putations. It seems to us that there is no way to quantify 
the error determined by possible deviations of the actual 
dielectric function from the Drude model. In our judg- 
ment, in order to obtain fully reliable predictions of the 
Casimir force with metallic plates it is necessary to find 
means of reducing the impact of the Drude extrapolation. 
In principle the straightforward way to achieve this goal 
would be by further extending optical data on the low 
frequency side. It is unlikely however that much progress 
can be made in this way, since it has been estimated 
that a determination of the Casimir force between two 
gold plates with an error of half per cent at a separation 
of 150 nm solely on the basis of measured optical data, 
would require extending optical measurements till wave- 
lengths in the millimeter region, which is a very difficult 
thing to do. 

An alternative approach, which does not require fur- 
ther extension of optical measurements beyond the limits 
of currently available optical data, was proposed recently 
by the author in Ref. [2^ . The idea is to introduce in the 
dispersion relations used to compute e(«^) suitable ana- 
lytical weights /(z), called by us "window" functions, 
that suppress the contribution of frequencies for which 
no optical data are available. We should observe that 
consideration of weighted dispersion relations is well doc- 
umented in the optical literature, where they have been 
exploited for example to derive new important relations 
connecting the refraction index to the extinction coeffi- 
cient. For details, we address the reader to Chapter 3 
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of ^2^. In our first work [25*1 we demonstrated that by 
taking weights f{z) that vanish both at zero and at in- 
finity, it is indeed possible to suppress very eff'ectively 
the contribution of frequencies outside the data interval 
Wmin < a; < Wmax- A simple family of window func- 
tions, parametrized by two integers p and q (with p < q) 
and by a complex frequency w, was offered there, and it 
was shown that for suitable choices of the parameters the 
error made by truncating the integral to the frequency 
range LUmin < < Wmax can be made negligible at both 
ends of the integration domain, in such a way that a pre- 
cise determination of e(j^) is possible, without extrapo- 
lating the optical data outside the interval where they are 
available. The key feature of the generalized dispersion 
relation that makes this possible is that, differently from 
the standard KK relation which only involves e"(w), the 
generalized relations involve both e'{uj) and e"{uj)^ and 
therefore they exploit the full information delivered by 
optical data. This is not a problem in principle if one 
considers that optical techniques exist, like ellipsometry, 
that permit to measure independently both e'(a;) and 
e"(w) with good precision. Such a technique was used 
for example in [22|. 

The analysis in [2^ was unrealistic however, in that 
it assumed that optical data in the interval Wmin < < 
Wrnax wcrc known with an infinite precision. Real optical 
data of course are not like that, as they are affected by 
ex per imental errors of various nature. That the analysis 
in [ 25l| was oversimplified has been pointed out recently in 
|15| . where the windowed dispersion relations were ten- 
tatively applied to the tabulated data for gold of Rcf. 
[20j . Using the same window functions and the same win- 
dow parameters used in , the authors of obtained 
negative values for the quantity e(j^) in certain intervals 
of the imaginary axis, and this is absolutely unaccept- 
able because on general grounds one knows that along 
the imaginary axis the electric permittivity of a causal 
medium must be positive (26j . In an attempt at under- 
standing their findings, the authors of Ref. 15j observed 
that the handbook [201 combines data from several dis- 
tinct experiments, using gold films prepared by different 
procedures. Since, as we said earlier, optical properties of 
gold films depend significantly on the deposition method, 
the operation of combining optical data for films having 
exceedingly different properties (that this is indeed the 
case with the tabulated data is shown in Sec. V below) 
may in principle spoil the overall KK consistency of the 
data, which is an essential precondition for dispersion re- 
lations to work well. We agree with the authors of Ref. 
p^ . though, that this potential problem with the tabu- 
lated data of Ref.fl^ cannot explain the obtained neg- 
ative values of ^{iC)- They argue that the main reason 
is the possible strong sensitivity of the windowed disper- 
sion relations to experimental errors in the optical data, 
in particular in the real part of the permittivity. If an 
exceedingly large amplification of small uncertainties in 
the optical data were an unavoidable feature of our win- 
dowed dispersion relations, their practical utility would 



be severely diminished, of course. The issue of error prop- 
agation is an important point deserving a detailed anal- 
ysis, and we address it in the present paper. We have 
performed a Monte Carlo simulation to determine how 
random errors in the optical data are propagated by the 
windowed dispersion relations. The results of these sim- 
ulation confirm the suspicion of the authors of Ref. [l^ , 
and explain their findings, revealing that the choice of the 
window functions that was made in [25| , was very unfor- 
tunate indeed, because it leads to a huge amplification of 
errors in the optical data. Trying to resolve the problem, 
we realized that the key property protecting the standard 
KK relation from this sort of instability is the fact that 
it involves a positive definite kernel. This property was 
badly violated by the weighted dispersions relations con- 
sidered in 25] , which are instead characterized by kernels 
attaining large negative and positive values. After sev- 
eral attempts, we could eventually fix this problem and 
we found a new class of weight functions, presented in 
this paper, that lead to positive kernels when applied to 
ohmic conductors. In this paper we demonstrate that 
the new weighted dispersion formulae do not suffer from 
large error propagation, and therefore they can be used 
to reliably compute e(i^), and in turn the Casimir force, 
solely on the basis of optical data extending from the 
IR to the UV, with no need of further extending optical 
measurements to lower frequencies. 

The plan of the paper is as follows: in Sec. II we 
review the KK formula routinely used to compute the 
dielectric function e(i^), and we discuss its weaknesses in 
the case of ohmic conductors. In Section III, we review 
briefly the weighted or "windowed" dispersion relations 
that were introduced in [25], and we present a new class 
of weights or " window functions" leading to dispersion re- 
lations with positive kernels when applied to conductors. 
In Sec. IV we present our numerical computations which 
demonstrate how the window functions considered in our 
first work suffer from a strong instability, and we show 
that the new weights resolve this problem. In Sec. V 
we apply the windowed dispersion relations to tabulated 
optical data for gold of Ref. ^], and in Sec. VI to some 
of the data in [22] . In Sec. VII we use the weighed dis- 
persion relations to derive an alternative expression for 
the so-called generalized plasma model that has been re- 
cently advocated as providing the correct description for 
ohmic conductors in the Casimir effect. Finally, Sec VIII 
contains our conclusions and a discussion of the results. 



II. DIELECTRIC FUNCTION AT IMAGINARY 
FREQUENCIES 

As it is well known from the theory of dispersion forces 
1] the Casimir interaction between two macroscopic bod- 
ies at temperature T depends on the dielectric functions 
of the two bodies, evaluated for a discrete set of imagi- 
nary Matsubara frequencies = 2TmkT/h, where n is 
any non-negative integer. Finding means of accurately 
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computing the quantities e(i^n) embodying the material 
dependence of the Casimir force, is the primary goal of 
this work. We shall see below that this is not a trivial 
task, in particular when ohmic conductors are consid- 
ered. As we said earlier, it is impossible to directly mea- 
sure e(i^) because any conceivable optical experiment can 
only determine the dielectric permittivity e{uj) along the 
real frequency axis. Knowledge of optical data along the 
real axis can be exploited to compute e(i^) by different 
mathematical procedures, depending on the sought accu- 
racy. For moderate precisions, at five per cent level or so, 
the simplest method is to make oscillator models of the 
optical data e{uj) for the material of interest. The explicit 
analyticity of such analytical models then permits to ob- 
tain the quantities e(«^Ti) by direct substitution uj — >■ 
into the oscillator formula. However, this straightforward 
procedure is not accurate enough for the purpose of inter- 
preting the most recent precision Casimir experiments. 
When higher precision is needed, the quantities e{i£,n) 
must be determined directly on the basis of optical data, 
without making recourse to simplified analytical models 
for the latter. This can be done by exploiting disper- 
sion relations entailed by causality of the electric permit- 
tivity, which relate the unmeasurable quantities e(i^n) 
to the measurable real- frequency permittivity e{uj). The 
dispersion relation that has been utilized until now in the 
literature is the following KK relation 



uje"{uj) 



(1) 



which is valid for all materials, like insulators and ohmic 
conductors, whose electric permittivity diverges at most 
like 1/uj at zero frequency. 

As we explained in the Introduction application of the 
above formula to ohmic conductors meets with difficul- 
ties. We recall first that the important range of Mat- 
subara frequencies that need be considered in Casimir 
computations depends on the separation a between the 
two bodies. In practice [2j], for an accurate determination 
of the Casimir force it is sufficient to consider Matsubara 
frequencies up to a maximum value of about ten times 
the characteristic frequency Wc = 2/ (2a) . For separa- 
tions a larger than 50 nm, which are the ones probed by 
current experiments, is is therefore sufficient to consider 
Matsubara modes with frequencies ^ ranging from 0.16 
eV/h, representing the frequency of the first Matsubara 
mode at room temperature, up to a maximum frequency 
of ten eV/h or so. Now we see from Eq. ^ that in 
order to evaluate e(z^) at any imaginary frequency ^ it 
is in principle necessary to know e"(a;) at all frequencies 
Lj. Unfortunately such a complete knowledge of e"(w) 
is never possible, because optical data are always re- 
stricted to some finite frequency range Wmin < uj < Wmax, 
starting from a non-zero minimum frequency Wmin > 0. 
Consider for example the data for gold collected in the 
handbook !20], which have been routinely utilized to in- 
terpret Casimir experiments using gold plates, including 
the recent experiments 11]. These data, which we shall 



describe more accurately later on, cover a range of fre- 
quencies extending from 0.125 eY/h, corresponding to 
the IR, till 10* cV/Ti. Unfortunately, this range is not 
sufficiently wide to permit an accurate estimate of the 
integral on the r.h.s. of Eq. ([T]) in the relevant range of 
imaginary frequencies ^, which as we said extends from 
0.16 eV/h to ten eV/h. In fact there is no difficulty 
on the high frequency side, because fall-off properties of 
e"{uj) entail that for all ^'s in this range, real frequencies 
UJ larger than a few tens of eV/h give already a negligi- 
ble contribution to the integral on the r.h.s. of Eq. ([T]). 
On the low-frequency side, however, we have a real prob- 
lem originating from the 1 /uj singularity displayed by the 
imaginary part of the permittivity of ohmic conductors. 
As a result e"(w) becomes extremely large at low frequen- 
cies, in such a way that the integral in Eq. ([1} receives 
a very large contribution from low frequencies for which 
no optical data are available. Since truncation of the in- 
tegral to the frequency Wmin results in a large error, one 
is forced to extrapolate the dielectric function e"(w) to 
frequencies uj < Wmin, where optical data are not avail- 
able, to evaluate the integral for uj < Wmm- As a rule, 
the extrapolation is done using the simple Drude model 



eDr(w) = 1 



uj{uj + ij) 



(2) 



where ujp is the plasma frequency, and 7 is the relax- 
ation frequency. Typical values for these parameters for 
gold are = 9 eV/h and 7 = 35 meV/?i [l^. As it 
is well known, the Drude model is expected to provide 
a reasonable approximation to the permittivity of ohmic 
conductors for low frequencies, and the hope is that the 
error resulting from use of the Drude extrapolation is not 
too large. According to this procedure, the quantity e(iO 
in Eq. ([l} is computed as 



e(iC) = 1 + ecut(«C) + eoxpt(iO , 



(3) 



where ecut (iS.) is calculated using the Drude extrapolation 
for e(uj) in the unaccessible frequency range uj < Wmin, 
while ecxpt(iC) is calculated using the experimental opti- 
cal data eg'xpt(w), according to the following formulas: 



2 p"" Loe'{,^iuj) 

Ecut «0 = - / duj 

Jo ^ + ? 



and 



^GXpt 



2 f 

= - / 



a;2 + ^2 



(4) 



(5) 



Later on, we shall denote by a subscript KK estimates 
for the permittivities and Casimir forces obtained by us- 
ing Eqs. (jSHSJ. The large weight of the Drude term in 
Eq. ^ has been clearly recognized recently [2^. Using 
the handbook data, it has been estimated (see Fig. 10 of 
[23 | and Fig. 2] below) that in the interval of imaginary 
frequencies ^ from 0.1 to 10 eV/U the quantity ecut(*C) 
gives the dominant contribution to e(i^) for ^ < 4 eV/h, 
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representing over 90 % of the total magnitude of e{i^) 
for ^ = 0.1 eV/Ti. We remark that the situation does 
not improve too much if the data interval is extended 
to longer wavelengths. For example, the authors of Ref. 
(23 | measured by ellipsometry the dielectric functions of 
several gold films for wavelengths from 33 to 0.14 /im, 
corresponding to the frequency interval from about 35 
meV/?l to 10 meV/h. Despite the fact that IR measure- 
ments were extended till photon energies that are roughly 
one fourth smaller than the handbook data, still ecut(*C) 
provides a large fraction of e(«C)i about 60 % for ^ = 0.1 
eV/fi and about 40 % for ^ = 1 eV/h (see thick solid line 
of Fig. H. 

To give a further sense of the importance of the Drude 
extrapolation for the Casimir effect, we have estimated 
the error which is made in the Casimir pressure P{a, T) 
between two gold plates if one neglects altogether the 
contribution of the Drude extrapolation in Eq. As it 
is well known, the Casimir force per unit area -P(a, T) 
between two identical plane-parallel homogeneous and 
isotropic plates, placed in vacuum at a distance a is pro- 
vided by the following Lifshitz formula [131 (a minus sign 
corresponds to an attraction between the plates): 



3or7\ 



p(a,r) = -M:^'^ 

n>0 a 



dk±k±q„ 



(6) 

where the prime over the n-sum means that the n ~ 
term has to be taken with a weight one half, a = TE, TM 
is the polarization, k± denotes the magnitude of the pro- 
jection of the wave- vector onto the plane of the plates and 
Qn = •y/fc^+^lT'^i where ^„ — 2TTn ksT/h are the Mat- 
subara frequencies. The quantities ra{iS,n, k±) denote 
the familiar Fresnel reflection coefficients of the slabs for 
Qf-polarization, evaluated at imaginary frequencies z^„. 
They have the following expressions; 



rTEii^n, k±) 



^(^Cn) Qn kn 

e(«6i) <?« + kn ' 



(7) 



(8) 



where fc„ — \/k\ -|- e(»C")?n/c^- Unless explicitly stated 
otherwise, it is understood that all Casimir pressures here 
and below are computed using the so-called Drude pre- 
scription [2j for ohmic conductors. In this scheme, the 
n = mode for TE polarization contributes nothing, and 
the ri > terms are evaluated using the dielectric func- 
tion e(«^) of the considered metal, without any modifica- 
tion aimed at suppressing the effect of ohmic dissipation 
for conduction electrons. In Fig. [l]we plot the fractional 
difference, in per cent, between the magnitudes |P(a,T)| 
and iPcutCa, T)| of the Casimir pressure between two par- 
allel gold plates at room temperature, that are obtained 
respectively if in e(i^) one includes or neglects the con- 
tribution of the Drude extrapolation ecut(*C) Eq. ([3]). 
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separation (^im) 



FIG. 1: Per cent error in the theoretical magnitudes of the 
Casimir pressure between two gold parallel plates at room 
temperature, resulting from neglect of the low-frequency con- 
tribution £cut(*0 i'l Eq. ((31 . The dashed line was computed 
using the handbook data of Ref.^^l ('^min = 0.125 eV/fi)- 
The solid line uses the optical data for sample 5 of Ref. [2^] 
(a;min = 0.042 eV/ft). 



The dotted curve displayed in the Figure is relative to the 
handbook data [205, for which Wmin = 0.125 eV /fi. Ex- 
trapolation to low frequencies was done using the com- 
'monly used values for the Drude parameters quoted in 
[l^l, namely — Q eV/h and 7 = 35 meV/?i. The solid 
curve in Fig. [T] was instead computed using the optical 
data relative to sample 5 of Ref. 



22|, consisting of an 



annealed gold film deposited on a mica substrate, with 
a thickness of 120 nm. We have selected this sample, 
because of the five samples studied in Ref. \22] this one 
has the closest plasma frequency {ujp = 8.38 eV/h) to 
the commonly accepted value for gold. In (2^ the op- 
tical data for this film were measured in the frequency 
interval extending from Wmin = 0.042 eY/h to Wmax = 9 
eV/U. The Drude parameters for this film were estimated 



in [231 to have the values ujr, 



.38 ± 0.08) eY/h and 



7 = (37.1 ± 1.9) meV/h. When considering this film, the 
handbook data from [20| were used to extend the data of 
[2^ above 9 eV/h. We see from Fig. [T]that at all sep- 
arations the Drude extrapolation has a smaller weight 
when the data of [l^l are used. This is a result of the 
fact, noticed already, that optical data of [l^] extend to 
lower frequencies than the tabulated data. We see from 
the Figure that the influence of the Drude extrapolation 
on the predicted pressure is very significant for both the 
handbook data and for the data of [2^, especially at 
submicron separations. For a separation a of 100 nm it 
contributes about thirty per cent of the total pressure in 
the case of the handbook data, and more than fifteen per 
cent for the data of Ref. [2J] . 

In the recent literature [ill, it has been claimed 
that relying on the handbook data of [l^l and by follow- 
ing the procedure outlined above, the magnitude of the 
Casimir pressure between two parallel gold plates can 
be predicted theoretically with an uncertainty of 0.5 % 
in the separation range from 150 nm to 400 nm. This 
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claim does not seem warranted to us, because of the large 
weight of the Drude extrapolation in the computation of 
the Casimir force. Considering the magnitude of this 
contribution, a precision of half per cent at separations 
around 200 nm can be claimed only if one is confident 
that the Drude extrapolation can be used to estimate the 
contribution to the Casimir pressure of unaccessible low 
frequencies w < Wmin with a fractional precision better 
than one or two per cent, a non trivial statement indeed. 
In our opinion, the reliability of Casimir computations 
would greatly increase if we could find means of substan- 
tially reducing the weight of the Drude extrapolation. We 
discussed earlier that one cannot expect to achieve this 
goal by further extending to lower frequencies the data 
interval. We shall demonstrate below that this goal can 
be achieved instead by using alternative forms of disper- 
sion relations, which strongly suppress the weight of the 
Drude extrapolation in determining the dielectric func- 
tion e(i^). 



fast in the origin and at infinity, it is in principle possi- 
ble to strongly suppress the contribution to the integral 
on the r.h.s. of Eq (|9]) of frequencies outside the inter- 
val [i^minjWmax] for which no Optical data are available, 
in such a way that an accurate estimate of £{i^) can be 
obtained by simply truncating the integral on the r.h.s. 
of Eq. ^ to the interval [wmim'^max]- We observe that 
by changing the form of the window function f{z) in Eq. 
^ , we have a means of changing at will the weight of the 
various spectral regions, and we can modify as well the 
relative weights of the real and imaginary parts of the di- 
electric permittivity. It is interesting to observe also that 
the windowed relations offer a possibility of checking the 
quality of the optical data, and in particular their degree 
of KK consistency, by verifying whether the obtained val- 
ues of e(iO do not change when the window function is 
changed. 

A. Old form of the window functions 



III. WEIGHTED KRAMERS-KRONIG 
RELATIONS 

In the work {25| we demonstrated that in principle the 
influence of the Drude extrapolation can be much de- 
creased by considering alternative dispersion relations to 
Eq. ([T|) involving appropriate weight functions, which are 
aimed at suppressing the contribution of low frequencies 
uj < Wmin, for which optical data are not available. 

Let us see briefly how the more general dispersion rela- 
tions come about. One starts from the following relation, 
generalizing Eq. ([T]), that can be easily proved by the 
contour integration method: 

(9) 

The above relation holds for any weight function /(z), 
later referred to as "window" function, that is analytic 
in the upper complex plane, and which satisfies there the 
following symmetry property 

f{-z*) = riz) , (10) 

In addition, the function /(z) and the permittivity e(z) 
should be such that the quantity u{z) = f{z){e{z) — 1) 
has at most a simple pole in the origin, and vanishes at 
infinity faster than 1/z" for some a > 0. The standard 
Kramers-Kroning relation Eq. ([T]) is a special case of Eq. 
(UJ, corresponding to the choice f{z) = 1. The distinc- 
tive feature of the general dispersion formula in Eq. ([9|l. 
as contrasted to the standard KK expression Eq. ([T]), is 
that in general it involves both the real and the imaginary 
parts of the permittivity. This is not a problem, in prin- 
ciple, because both quantities can be measured, using for 
example ellipsometry. The key remark made in (25l | was 
that by using window functions that vanish sufficiently 



As an example, in [25| we considered "window" func- 
tions f{uj) of the following form 

= (11) 

where w is an arbitrary complex frequency such that 
Im(z/;) < 0, and p and q are non-negative integers such 
that < p < q. A practical application of the above win- 
dow functions was made in Ref. [l^, where they were 
used to compute the permittivity e(i^) of gold, on the 
basis of the tabulated handbook data [1^. The used 
windowed parameters were the same that had been used 
earlier in ^2^ i.e. w = (1— 2j) eV/h, p — 1 and q = 3. The 
obtained results were very unsatisfying, since large devi- 
ations from the KK results were found and, worse than 
this, unacceptable negative values for e(i^) were found for 
values of ^ in the interval from 2.44 to 2.92 eV/h, and 
again for ^ > 7.8 eV/h. The authors of fis'l commented 
that the problem was presumably determined by an ex- 
ceeding propagation of experimental errors affecting the 
optical data by the windowed relations. 

After closer inspection we realized that the problem 
originated from the inconvenient choice of the analytic 
form of the window function Eq. (Ilip that we made in 
our first work. While a more systematic investigation 
of the problem is postponed until the next Section, we 
can easily understand what is wrong with the window 
functions in Eq. (jlip by considering the quantity g{uj) 

- ^^-^iHfi^mu:) - 1)] , (12) 

whose integral along the positive real axis should repro- 
duce, according to Eq. ([9]), the permittivity e(«^)-l. In 
Fig. [2] we plot the function g{uj), using for e(w) the hand- 
book data. The displayed curve corresponds to ^ = 2.5 
eV/?i, and to the window parameters p = 1, g = 3 and 



7 



1000 

500 


-500 
-1000 
-1500 
-2000 




4 6 
C(eV/fi) 



10 



FIG. 2: Plot of the quantity g(ijj) relative to the handbook 
optical data for gold of Ref. for ^ = 2.5 eV/ft. The 

window function f{z) is as in Eq. (jlip . for window parameters 
p = 1, g = 3 and w — {1 ~ 2i) eV /h. 



w = (1 — 2i) eY/h. The chosen value of ^ is one of those 
for which negative values of e(i^) were found in Ts*]. We 
see from the Figure that g{uj) attains large positive and 
negative values. This implies that even slight errors in 
the optical data may strongly affect the delicate balance 
between positive and negative regions in the integral on 
the r.h.s. of Eq. leading to large errors in the ob- 
tained value of e{i^). That the estimates of e(i^) are 
indeed very sensitive to errors in the optical data will 
be demonstrated more thoroughly in the next Section. 
Here, we remark that this instability can partially cured 
by using different values for the window parameter w, 
such that the quantity ^(w) displays lower peaks. How- 
ever, we prefer to look for a more radical solution of the 
problem, by considering alternative forms of the window 
functions, which are presented in the next Section. 



B. Improved choice of the window function 

A highly desirable feature displayed by the standard 
KK relation Eq. ([T|) is that the integrand on its r.h.s. 
is a positive definite quantity, since the imaginary part 
e"(w) of the permittivity of all materials is positive. This 
property of the integrand is essential in ensuring robust- 
ness of the integral with respect to small errors in the 
optical data, which permits to estimate reliably the per- 
mittivity e(i^) along the imaginary frequency axis, pro- 
vided only that data are available in a sufficiently large 
frequency interval. This led us to wonder whether there 
exist choices of the window function f(z) such that the 
integrand on the r.h.s. of the generalized relations Eq. 
^ has a definite sign, in such a way that the estimate 
of becomes more robust, while at the same time 

preserving the key feature of the generalized dispersion 
relations of suppressing the contribution of low frequen- 
cies, for which no optical data are available. Achieving 
this goal is not so easy because differently from the stan- 
dard KK relation, the integrand on the r.h.s. of Eq. © 



has no definite sign in general, as it can be seen for ex- 
ample in Fig. [21 This is so because for window functions 
f{z) different from one, the integrand in Eq. © involves 
both the real and the imaginary parts of the permittivity 
e(a;), and it is well known that differently from e"(w), 
the real part of the electric permittivity e'(w) has no def- 
inite sign. In addition to this, the quest for an integrand 
of definite sign is further complicated by the fact that 
analyticity prevents the real and the imaginary parts of 
the window functions from having a definite sign, if one 
insists that they should vanish both at the origin and at 
infinity (see discussion in Ref. [2^ ) . 

Notwithstanding these general complications, fortu- 
nately enough we managed to find a class of window 
functions for which the integrand on the r.h.s. of Eq. 
([5]) is in fact positive, when ohmic conductors are con- 
sidered. In order for this to be possible, we had to relax 
the condition that the window functions vanish at infin- 
ity, and instead consider weight functions that approach 
one at infinity. This is not a problem at all for practical 
purposes, because as we said earlier optical data extend 
to sufficiently high frequencies for no further suppression 
to be necessary with the help of window functions on the 
high frequency side. The window functions that fit our 
needs have the following simple expression: 



f{z;b) 



- 62 



(13) 



where b is an arbitrary non- negative real frequency 6 > 0. 
As we see these functions vanish in the origin, which en- 
sures the desired suppression of low frequencies in the 
integral on the r.h.s. of Eq. However, differently 

from the previous window functions in Eq. (|lll) . these 
functions enjoy also the nice feature that their real and 
imaginary parts are both semi-definite functions along 
the positive frequency axis, since their imaginary part is 
negative for lo < b and vanishes for uj > b, while their 
real part is zero for ui < b and positive for uj > b. In- 
terestingly, weighting functions of the form (cj^ — 52^-1/2 
have been used in the past to obtain relations connecting 
the refraction index n{uj) in the interval [0, b] to the ex- 
tinction coefficient k{uj) in the remainder of the spectrum 
(see Chapter 3 in [20|). Insertion of Eq. ([T5)) into Eq. 
© results into the following formula for e(z^): 



.(.o = i+|Ji+Q 



did- 



l-e'{u:) 



^2 + e 



e"{u) 



52 



(14) 



We observe that the formula involves only the real part of 
the permittivity e'(w) for w < 6, and only its imaginary 
part f" {uj) for uj > b. The integrand involving e"(w) on 
the r.h.s. of the above formula is positive definite, while 
the integrand involving e'{uj) is positive whenever e'(w) is 
less than one. This is what we needed, because the real 
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FIG. 3: Plot of e'M - 1 for gold. The solid line is for the 
handbook data of [201 ] . the dashed one for sample 5 of Ref. 

m. 



part of the conductivity of ohmic conductors is charac- 
terized by the fact of being less than one essentially at 
all frequencies, and surely below the interband transition 
'Winter, as it it Can be seen from Fig. [31 Therefore, if we 
choose b to be any frequency such that e'(a;) is less than 
one for uj < b both integrands in Eq. ([14p are positive, 
similarly to the standard KK relations. Contrasted with 
the standard KK relations, the new relation Eq. ([Til) re- 
ceives a much smaller contribution from low frequencies. 
This is so because in the limit of vanishing frequency 
the integrand in the first term of Eq. dTTl) vanishes like 
w^, since e'(w) stays finite in the static limit, while the 
integrand of the standard KK relation approaches a non- 
vanishing constant in the same limit, due to 1/uj singu- 
larity displayed by e"(w). Thus we may expect that the 
new relation permits to estimate e(«^) more reliably than 
the standard KK relation because, while enjoying the ro- 
bustness of a positive kernel, the weight of the Drude 
extrapolation is much smaller. In practical applications 
of Eq. ([T4l) we shall typically consider frequencies b that 
belong to the so-called non-relaxation region 6^7, and 
are smaller than the interband transition Winter- For such 
values of b, the conditions e'(a;) < and le'(iw)l ^ 1 hold, 
and therefore the values of e'{uj) should not be affected by 
large experimental errors, contrary to what may happen 
at frequencies for which e'(w) nearly vanishes ]15]. 

For the purpose of numerical evaluation, the occur- 
rence of an integrable singularity for lo — b in the inte- 
grals on the r.h.s. of Eq. (ITTl) is quite inconvenient. We 
can easily dispose of this singularity, however, by per- 
forming the change of variables uj = b sin y in the first 
integral on the r.h.s. of Eq. ([TTl) . and uj — b coshy in the 
second one. After doing this, Eq. transforms into: 



em 




7r/2 

dy 



sin^ y 



: [l-e'(6siny)]+ / dy 
Jo 



cosh^ y 



cosh^ y I 



'(6 coshy) 



(15) 



IV. NUMERICAL COMPUTATIONS 



We have tested the performance of the windowed dis- 
persion relation using the improved choice of the weight 
function Eq. ([T4l) on some of the data sets for gold that 
have been used recently for Casimir computations. First 
we consider the handbook data [1^. Since they have 
been widely used to interpret Casimir experiments, in- 
cluding in__particular the short-separation precise exper- 
we shall reserve special attention to these 
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iments 

data, and we shall spend time to evidence some incon- 
veniences of these data that must be handled with some 
degree of caution. In addition to the handbook data, we 
used the data of the recent work [l^j. In principle, these 
data present several advantages over the handbook data. 
First of all, they were obtained by using ellipsometry, 
which is a powerful optical technique which permits to 
measure independently the real and the imaginary parts 
of the dielectric function. Since our dispersion relations 
explicitly involve the real part of the permittivity, it is 
important to have good quality data for this quantity. 
In addition to this, the data of [22| present two desirable 
features: on one hand, as we observed earlier, the data in 
[2^ extend down to a frequency a;min=38 meV/h, which 
is lower than the minimum frequency a;min=0.125 eV/h 
reached by the handbook data. We shall find indeed 
that a lower frequency around 40 meY/h is sufficient to 
compute with high precision the Casimir force, using the 
windowed dispersion relations. The second and perhaps 
more important feature of Ref. [22I] is that it provides 
homogeneous data on well defined gold samples over the 
wide range of wavelengths from 0.14 to 33 /um, which give 
the main contribution to the Casimir effect. This is not 
so for the handbook [g^l, whose data in the important 
spectral region from the near IR to the near UV combine 
results from two different experiments [29l Isoj , utilizing 
differently prepared gold films. We shall see later on that 
the optical properties of these gold films display striking 
differences, and therefore wrong results may easily result 
if the handbook data are used in computations where KK 
consistency of the data plays an important role. 



A. Suppression of the Drude contribution 

Our key motivation for introducing windowed disper- 
sion relations is that they strongly reduce the weight of 
the Drude extrapolation, and correspondingly increase 
the weight of experimental optical data, in such a way 
that obtained values for e(i^) are much less sensitive to 
experimental uncertainties in the Drude parameters. We 
can verify to what extent the new windowed relation Eq. 
(ITTl) suppresses the contribution of the Drude extrapo- 
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FIG. 4: Relative contribution of the Drude extrapolation 
ecut(iO to the dielectric function according to the stan- 

dard KK formula (thick lines) and to the windowed formula 
Eq. (|14p (thin lines). Solid lines were computed using optical 
data of sample 5 of Ref. f2^ (solid lines), while the dashed 
lines were computed using the handbook data [20|] . 

lation by proceeding in a way similar to Ref. We 
decompose the quantity e(i^) in Eq. (121) in a way analo- 
gous to Eq. (O: 

eiiO = 1 + ecut(i^; b) + eexpt(i^; b) , (16) 

where ecut(*C!^) is calculated using the Drude extrapo- 
lation for e(w) in the unaccessible frequency range uj < 
Wijiijj, while ecxpt(*Ci b) is calculated using the experimen- 
tal optical data, according to the following formulas: 

<^cut{i£,;b) = 

^7(^ I d-^:^Mn-; b){eM-) 1)] , 

(17) 

and 

eoxpt(j^;6) = 

—frrJ dc.^--2lm[/(c.;6)(eexptH-l)] , (18) 

where f{z\ b) is as in Eq. (IT^ . In order to demonstrate 
how well the windowed formula suppresses the Drude 
contribution ecut(*Oj Fig- SI we plot the relative con- 
tribution of the Drude extrapolation f-cut{iC) / ^i^C) ^ that 
obtains if one uses the standard KK relation (thick upper 
solid and dashed lines) or the generalized dispersion rela- 
tion (thin lower solid and dashed lines) with f{z;b) given 
by Eq. (fT3| . for 6 = 1 eV/h. The displayed imaginary- 
frequency range extending from 0.1 eY/h to ten eY/h 
is the one relevant for determining the Casimir force for 
plates separations larger than 100 nm. The solid lines in 
Fig. 0] were computed using the optical data relative to 
sample 5 of Ref. [22| , while the dashed lines are relative 
to the handbook data HO]- We recall that in ^23] optical 



data were measured in the frequency interval extending 
from Wmin — 0.042 eY/U to Wmax = 9 eY/h, and so for 
UJ larger than 9 eY /fi we used data from Ref. [l^l to 
compute ecxpt(iO- The Drude parameters used to com- 
pute ecut(*0 were as follows: for sample 5 of j2^] we used 
the average values quoted there, i.e. oOp = 8.38 eY/h 
and 7 = 37.1 vcieY /h, while for the handbook data we 
used the commonly used values ojp = 9 eY/h and 7 = 35 
meV/?i. It is apparent from Fig. |4] that the general- 
ized dispersion relation is very effective in suppressing 
the contribution of the Drude extrapolation, in compar- 
ison to the standard KK relation. This is especially so 
for the data of (22)], due to the fact that they extend to 
smaller frequencies than the handbook data. In this case, 
for the lowest displayed imaginary frequency of 0.1 eY/h, 
the relative contribution of ecut is reduced from 70.2 % to 
8.4 %, but already for ^ = 0.5 eY/h it is reduced from 55 
% to 1.5 %, and it becomes rapidly negligible for larger 
imaginary frequencies. 

We can show that thanks to suppression of the Drude 
extrapolation, experimental uncertainties in the Drude 
parameters lead to much smaller uncertainties in the di- 
electric function e{i^), when the window method is used 
in comparison with the standard KK formula. This in 
turn implies that the Casimir force can be predicted more 
accurately by the widow approach. This is an important 
result because experimental uncertainties in the Drude 
parameters have a non-negligible impact on the Casimir 
force, when standard KK relations are used. To be defi- 
nite, we considered the optical data of [2^. In this work 
the Drude parameters for the used gold films were deter- 
mined by several methods, and precisely by fitting with 
the Drude formula both e(a;) and the complex refraction 
index n(cLi) = -y/ e(aj), and also by performing a KK anal- 
ysis of the data. Depending on the used method, the 
obtained values of ujp and 7 for the five available gold 
samples showed non-negligible variations, ranging from a 
minimum of 1 % to a maximum of 2 % for ujp, and from 
5 % to 14 % for 7. It was found that by itself this un- 
certainty in the Drude parameters entails an uncertainty 
in the Casimir force that can reach 1 % at separations 
around 100 nm. In order to demonstrate the effectiveness 
of the window method in reducing this source of uncer- 
tainty, in Fig. Owe plot the relative per cent variation of 
the dielectric function e{i^) determined by small changes 
in the Drude parameters, according to the standard KK 
formula (thick lines) and to the windowed formula Eq. 
([T4| (thin lines). The upper four lines correspond to a 
2 % increase in Wp, while the four lower lines correspond 
to a 14 % increase of 7. Solid lines were computed using 
optical data for sample 4 of Ref. 22] (solid lines) , while 
the dashed lines were computed using the handbook data 
pol ]. In Fig. [6] we display the corresponding per cent vari- 
ations in the Casimir pressure between two gold plates 
in vacuum, at room temperature. The lines in Fig. [6] 
have the same meanings an in Fig. [SJ apart from the fact 
that now the 4 lower curves correspond to a 2 % increase 
in ujp, while the four upper curves correspond to a 14 % 
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FIG. 5: Relative per cent variation of the dielectric function 
e(i^) determined by small changes in the Drude parameters, 
according to the standard KK formula (thick lines) and to the 
windowed formula Eq. (|14|l (thin lines). The upper four lines 
correspond to a 2 % increase in Up, while the four lower lines 
correspond to a 14 % increase of 7. Solid lines were computed 
using optical data of sample 4 of Ref. [2^ (solid lines), while 
the dashed lines were computed using the handbook data . 



FIG. 6: Relative per cent variation of the Casimir pressure 
between two gold plates at room temperature, determined by 
small changes in the Drude parameters, according to the stan- 
dard KK formula (thick lines) and to the windowed formula 
Eq. H14|) (thin lines). The lower four lines correspond to a 2 
% increase in LUp, while the four upper lines correspond to a 
14 % increase of 7. Solid lines were computed using optical 
data of sample 4 of Ref. j22l] (solid lines), while the dashed 
lines were computed using the handbook data 20j. 



increase of 7. We clearly see that uncertainties in the 
Drude parameters have a much smaller impact when the 
window approach is used, as compared to the traditional 
KK relations. It is also clear that the window method is 



more effective when applied to the data of [22| than for 
the handbook data. Indeed we see that when the former 
data are used, the maximum uncertainty in the Casimir 
pressures decreases from about 1 % to 0.05 %, while for 
the handbook data the maximum uncertainty decreases 
from 0.75 % to 0.14 %. Similarly to Fig. H the supe- 
rior performance of the window method when the data 
of [22J are used is a consequence of the fact that these 
data extend to smaller frequencies than the handbook 
data. 



B. Propagation of random errors 

The next thing that we tested is how possible errors 
in the optical data are propagated by windowed disper- 
sion relations. Indeed, if we could rely on perfect optical 
data with no uncertainties, the mathematical properties 
of analytic functions would permit to compute the ex- 
act values of the permittivity anywhere else in the upper 
complex plane, starting from knowledge of the dielectric 
permittivity in any however small interval of the real fre- 
quency axis. The existence of experimental errors in the 
optical data makes this goal impossible to achieve, and 
in practice if one tries to compute the dielectric function 
at points of the complex plane that are too far for the 
experimental data, one inevitably ends up with estimates 
having large errors. In view of this practical limitation, it 
is important to investigate quantitatively how uncertain- 
ties in the optical data are propagated by the dispersion 
formulae used to compute e(i^). This is not an easy ques- 



tion to address, because a complete answer can be given 
only if sufficient information is provided on the nature 
the experimental errors affecting the data. 

Since neither the handbook ^2d\ nor Ref. [l^] provide 
detailed information on the magnitude of the experimen- 
tal errors affecting their data, we performed a Monte 
Carlo simulation to estimate how statistical errors in the 
optical data are propagated by the windowed dispersion 
relations, under the reasonable assumption of a statisti- 
cal uncertainty of a few per cent in the optical data. We 
proceeded as follows. Optical data are usually presented 
as lists of values S = {fT.oxp(wi), fccxp('j-'i)}, i = 1,...N 
for the refraction index n{io) and the extinction coeffi- 
cient fc(cLi) {n{uj) + i k{uj) = -y/ e(aj)), for a discrete set of 
frequencies Wmin = uii < UJ2 < ■■■ < llin = i^max- 
For any definite set S of data, say the handbook data, 
we computed the quantity eexpt(iC) using the general- 
ized formula in Eq. (fT8|) . We then made the simple 
assumption that the experimental values of n and k in 
S had a common statistical per cent error Scxp- In order 
to estimate the corresponding error (5ecxpt(*0: '^^ ^'^'H- 
domly generated M new hypothetical sets of data Sa — 
{(n^"\ fcf ^), • • • , {nP,kP)}, a = 1, 2, • • • M. Each data 
set Sa was generated by extracting 2N random numbers 
{rii, ki} , i = 1 . . . A'^ from 2N independent gaussian dis- 
tributions, having mean values respectively equal to the 
experimental values {fT.oxp(i^i)j ^exp(wi)}, and variances 
respectively equal to {(ncxp(t^i) ^cxp)^, (fcexp(wi) <5exp)^}. 
Each data set Sa was then used to obtain a new estimate 
^cxpt(*C) s-nd the average (absolute) error 5eexpt(«0 ou 



11 



200- 



150- 



100- 



50- 




4 6 



10 



0.6- 



-H 

^ 



0) 



.5 



. 4 



0.3- 



4 6 
C(eV/ii) 



10 



FIG. 7: Simulated uncertainty (in per cent) in the value of 
the dielectric function e(i^) versus imaginary frequency (in 
eV/h). The dielectric permittivity e(i^) was computed us- 
ing the handbook data [S^l, with the help of the windowed 
dispersion relation for the choise of window functions in Eq. 
(|ll|l . The used window parameters are p — 1, q = 3 and 
w = {1 — 2i) eV /h. In the simulation we assumed a 3 % error 
in both n{uu) and k{ijj). 



FIG. 8; Simulated uncertainty (in per cent) in the value of 
the dielectric function e(i^) versus imaginary frequency (in 
eV/h). The dielectric permittivity e(i^) was computed using 
the handbook data [20l |. with the help of the windowed dis- 
persion relation Eq. (|14p . for b — 1 eV/h (solid line) and for 
b — 3 eV /h (dashed line). In the simulation we assumed a 3 
% error in both n{ij) and A:(aj). 



texpt 

(z^) was estimated by the formula: 



\ 



M 



M - 1 ^ 



■ (19) 



In our simulation we took M — 1000. We studied by this 
method both the old window functions Eq. (fTTj) , and the 
new ones Eq. ([T^. 

We recall that by using the windowed relations with 
window functions of the form given in Eq. (|11|) the au- 
thors of [3| obtained unacceptable negative values for 
e(i^) in certain imaginary frequency ranges, and also 
for the frequencies for which positive results were found, 
large deviations from the KK values were observed. In 
order to explain the findings of [isl] , we investigated first 
these window functions, for the same values of the win- 
dow parameters p — 1^ q ^ i and w = (1 — 2i)eV/h 
that were used in our earlier work [25|, and that were 
subsequently considered by Geyer et. al [l^. Results 
of the Monte Carlo simulations give support to the con- 
jecture made by the authors of [15|, that this choice of 
the window function determines an exceeding amplifica- 
tion of experimental errors in the optical data, explaining 
the obtained negative values for e(i^). This can be seen 
from Fig. [71 where we plot the fractional random error 
<^eoxpt(iC)/£KK(iC) (ill P6r cent) for the above choice of 
window parameters, assuming a 3 % error in both n{uj) 
and k{uj). As we see, the error 5ecxpt(*0 very large at 
all frequencies. The largest error occurs for values of ^ 
near the zero = 2.4 displayed by the window function 
along the imaginary axis. It is clear that estimates of 
e(i^) obtained with this choice of the window function 
do not have much meaning. We verified that better re- 
sults can be obtained for different values of the window 
parameter w. For example, for w = —5i eV/h the un- 



certainty (5eoxpt(*^) becomes less than 7 per cent, in the 
entire range of ^ from 0.1 and 10 eV/h, again assuming 
a 3 % error in both n(uj) and k(uj). 

Even though it is possible to sensibly reduce the errors 
by suitably tuning the value of w, in a ^-dependent way, 
we did not pursue further the old form of the window 
functions, because much better results can be obtained 
by using the new window functions given in Eq. (1131) , and 
the associated dispersion relation Eq. (ITi)) . With such a 
choice of the window function, the uncertainty in ecxpt(*C) 
decrease tremendously. This is shown in Fig. |8l where we 
plot the fractional uncertainty 5eexpt(*C)/£KK(^C) (in per 
cent), versus the imaginary frequency ^ (in eV/h), for the 
two values of the window parameter b — 1 eV/h (solid 
line) and 6 = 3 eV/h (dashed line), again assuming a 3 
% error in the handbook data for n{u!) and k{uj). We see 
that the uncertainty ^ecxpt(iC)/^KK(iC) is now less than 
0.7 % for all displayed imaginary frequencies. Thus our 
Monte Carlo simulation gives support to the expectation 
that the new window functions permit to obtain robust 
estimates of the dielectric function e(iC), thanks to the 
positive definite character of the corresponding integrand 
on the r.h.s. of Eq. In view of these findings, we 

shall not consider any longer the old window functions 
in the next Sections, and we shall exclusively rely on the 
new ones for our next computations. 



V. AN APPLICATION TO TABULATED DATA 
FOR GOLD 

In this Section we shall use the new window functions 
and the associated dispersion relation Eq. ([Tl)) to obtain 
an estimate of dielectric permittivity e(j^) of gold, on the 
basis of the handbook data [20[ . 

Before we do it, an important observation on the con- 
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TABLE I: Values of the complex permittivity for gold films 
from Ref. [2^ and Ref. [S^], that have been both included in 
the handbook IM. 



cj (in eV/h ) 


e{uj) (Ref. [29J) 


e{uj) (Ref. [30]) 


0.6 


-168.2 + i 23.3 


-170.7 + i 25.6 


0.7 


-125.0 + i 15.6 


-165.9 + i 18.7 


0.8 


-96.0 + i 11.0 


-131.9 + i 12.65 


0.9 


-76.7 + i 7.96 


-94.3 + i 8.9 



sistency of these data is in order. The handbook [20| 
quotes values for the refraction index n{uj) and the ex- 
tinction coefficient fc(cj) of gold in the frequency range 
from Wmin = 0.125 eV/?i up to Wmax — 9919eV/?i. An 
important feature of these data is that they collect to- 
gether results from different experiments, performed in 
different spectral regions, utilizing gold films prepared 
by different procedures. In particular, data in the in- 
frared region of the spectrum, with frequencies in the 
interval 0.125 eV/?i < < 0.98eV/?i, were taken from 
Ref. (2^, which used an evaporated gold film on a pol- 
ished glass substrate. Data in the frequency region 
0.6eV/?i < 1^ < 6eV/?i were taken from Ref. [30| . 
which used annealed films evaporated in ultrahigh vac- 
uum on fused silica substrates. Data in the interval 
6.199 eV/?i < a; < 26.38 eV/?i were taken from Ref. 
[3ll | , which used evaporated films onto polished glass sub- 
strates, in a conventional vacuum system. Finally, data 
in the interval 26eV/?i < w < 88eV/?i were taken from 



Ref. [32|, which used thin gold films evaporated onto 
substrates of collodion. For the purpose of Casimir com- 
putations, the important spectral regions are the IR one 
covered by Dold's et al. data [1^ , and the region extend- 
ing from the near IR to the near UV in Theye's data [soj . 
At a closer inspection, one realizes that these two sets of 
data are quite inconsistent with each other, probably be- 
cause of the different deposition procedures adopted in 
these measurements. This can be appreciated by com- 
paring the values, presented in Table 1, of the dielectric 
function e{uj) computed from the two sets of data, in the 
spectral region where they overlap i.e. from 0.6eV/?i to 
0.9 eV/h. We clearly see that there exist large differences 
between the two samples. 

A more striking measure of the inconsistency between 
the two data sets is obtained by estimating the respective 
plasma frequencies Wp. This is an important quantity 
to consider, because as we saw earlier its value has a 
large impact on the magnitude of the Casimir force. We 
have estimated ujp by the following procedure. It is well 
known ^33] that for frequencies uj in the non-relaxation 
region lj ^ ^ but smaller than the interband transition 
Winter (^intor — 2 e\f/fi for gold), the real part of the 
electric permittivity of metals can be well described by 
the formula: 



-50- 
-100- 
-150- 
-200- 



4, 



FIG. 9: Real part of the dielectric function of gold in the 
infrared and red regions, versus square of wavelength (with A 
in microns). Displayed data are from Ref. [S^l (squares) and 
[29( 1 (diamond), which are part of the handbook [S^]. Also 
displayed are the linear fits to the data (solid and dashed 
lines), according to Eq. (|20|l . The used fit parameters for 
the plasma frequency are uip = 9.19 eV /fi (solid line) and 
ujp = 7.1 eV/h (dashed line line). 



where A is the wavelength and eintcr is the contribution 
of interband transitions, which for frequencies sufficiently 
below Winter can be considered as a constant. From the 
above relation, we expect that e' should vary linearly with 
the A^, with a slope proportional to ujp. We determined 
ujp by fitting the data in the infrared and red regions 
of the spectrum. Since Theye's data [s^l extend to the 
interband region, we did not fit parameter 

for these data, but we computed its value by means of 
the dispersion formula: 



Winter 



'-r ^."H, (21) 



1 H" Winter 



(20) 



from which we obtained einter = 5.9. In the case Dold's 
et al. data f29], einter was used as a fit parameter. We 
note however that the values of einter have little impact on 
the obtained values of the plasma frequency uj}, because 
for the considered wavelengths |e'(a;)| ^ ejntor- The fits 
resulted in the following values of the plasma frequency 
for the two sets of data: ojp = 9.19 eV/h for Theye's data 
and ujp = 7.7 eV/h for Dold's et al. data. The data and 
the fitting straight lines are displayed in Fig. ^ from 
which we see that the linear fits are very good for both 
sets of data. The obtained value of ojp for Theye's data 
compares well with the value of 9 eV/h, that was obtained 
on theoretical grounds by Lambrecht and Reynaud [13], 
as well as with the value of 8.9 eV/h that was obtained 
by the authors of [ll| from resistivity measurements on 
their gold samples. The value of Wp that we obtained for 
Dold's et al data is consistent with the value of 7.5 eV/h 
that was obtained by a similar method in Ref. [IJ . In the 
latter work however no attempt was made to determine 
the plasma frequency for Theye's data separately. 

Since the value of the plasma frequency for Dold's data 
is so much smaller than that for Theye's data, we have de- 
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FIG. 10: Relative contribution of the Drude extrapolation 
ecut(*0 to the dielectric function e(i^) of gold, for the stan- 
dard KK relation (thick dotted line) using the complete hand- 
book data (oimin = 0.125 eV/h), and the windowed dispersion 
formula after exclusion of Dold's et al. data [2^ (cJmin ~ 0.6 
eV/h) for b = 1 eY/h (thin solid line) and 6 = 1.5 eV/h (thin 
dashed line). 
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FIG. 11: Per cent difference between the window and the 
KK estimates for the permittivity e(if) of gold, computed 
using the handbook data [S^]. The windowed estimate was 
computed using Eq. H13|) . after excluding Dold's data, for 
b = 1 eV/h (solid line) and b = 1.5 eV /h (dashed line). The 
Drude parameters used to extrapolate the data are LUp = 8.9 
eV/h and 7 = 35.7 meV/h. 



cided not to use Dold's data in our computation of the di- 
electric function along the imaginary frequency axis. The 
minimum frequency Wmin of the remaining data therefore 
coincides with the one Wmin — 0.6 eV/h for Theye's data. 
The Drude parameters used to extrapolate the data be- 
low 0.6 eV/% were same as in [ll|, i.e. cOp = 8.9 eV/U 
and — 35.7 uieV/h. The price we had to pay for ex- 
cluding Dols'd data is that the contribution of the Drude 
extrapolation returns to be significant, even if the win- 
dowed dispersion relation Eq. (fT4|) is used, as it can be 
seen by comparing the thin lines in Fig. [TU]with the dot- 
ted thick line in Fig. |4l However, thanks to the ability of 
the windowed relation to suppress low frequencies, we see 
that despite the large value of w = 0.6 eV/h, the weight 
of the Drude extrapolation in the windowed relation is 
significantly less than the one (thick dotted line in Fig. 
[TU| resulting from the full set of tabulated data, begin- 
ning from 0.125 eV/h, when the traditional KK relation 
is used. 

We compared the values of e{i£,) obtained from Eq. 
(IT4l) with those obtained from the standard KK rela- 
tion Eq. ([T]). We shall denote the two estimates by 
ewin(*C) s-nd eKK(i'C)j respectively. In Fig. [Tl]wc display 
the fractional difference (ewin(*0 ~ eKK(*C))/^KK(i'C) (in 
percent) between the two estimates. The solid and the 
dashed lines are for 6 = 1 eY/h, and for b — 1.5 eY/h, 
respectively. We remark that the standard KK values of 
£kk(*C) in Fig- HH have been computed using the full 
handbook's data, and therefore they coincide with the 
values used in [lli| for comparison of the Drude model 
approach with experiments. We note that the windowed 
values of e(«^) obtained for the two chosen values of the 
window parameter b differ from each other by less than 
one per cent, and both differ from the KK values by 
roughly the same amount. So we see that the improved 
window method produces results that are substantially 
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FIG. 12: Per cent difference between the window and the KK 
estimates for the permittivity e(if) of gold, both computed 
using the full handbook data [2y]. The windowed estimate 
was computed using Eq. (|13p . with 6=1 eV/h (solid line) 
and b — 1.5 eV/h (dashed line). The Drude parameters used 
to extrapolate the data are ujp = 8.9 eV/h and 7 — 35.7 
meV/?i. 



consistent with the traditional method. We think how- 
ever that the window estimates obtained here are never- 
theless more reliable than the traditional ones, because 
they rely to a less extent on the Drude extrapolation. 

It is interesting to see how the above results change 
if we use the window method with the full set of hand- 
book's data. This can be seen from Fig. [T^] where we 
plot the corresponding fractional difference (ewin(iC) — 
eKK(*C))/^KK(iC) (in per cent), the solid and the dashed 
lines have the same meaning as in Fig. 1111 It is evident 
that after we include Dold's data the agreement between 
the window and the KK estimates worsens significantly. 
We interpret this fact as a further proof of the KK in- 
consistency that results from combination of Theye's and 
Dold's data. 
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FIG. 13: Per cent difference between tfie window and tlie KK 
estimates for the Casimir pressure between two gold plates, 
computed using the handbook data [20l ] . The windowed pres- 
sures Pwin were computed excluding Dold's data, by using 
Eq. ([131) for & = 1 eV/^i (solid line) and b = 1.5 eV/Ti (dashed 
line) . The Drude parameters used to extrapolate the data are 
ujp = 8.9 eV/h and 7 = 35.7 meV/h. 



We used the windowed estimates of e{i^) to compute 
the Casimir pressure between two parallel gold plates in 
vacuum. We let Pwin (a, 7") and PKK{a,T), respectively, 
the Casimir pressure obtained by plugging into Eq. ^ 
the windowed and the traditional KK estimates for the 
permittivity In Fig. [13] we plot the per cent dif- 

ference [Pwin (a, T) — PkkCo, T)]/PKK(a, T) between the 
windowed and the KK pressures computed using the 
handbook data, versus plate separation (in microns) . As 
before, the windowed values were computed excluding 
Dold's data, while the KK pressures use the full set of 
data. Two different values of the window parameter b 
were considered, namely b — 1 eV/Ti (solid line) and 
b — 1.5 eV /h (dashed line). As we see, the pressures 
computed by the window method differ from those ob- 
tained by the traditional KK method by less than 0.2 % 
in the entire displayed separation range. The width of 
the region comprised between the solid and the dashed 
lines in Fig. [11] provides an indication of the uncertainty 
in the windowed prediction of the Casimir pressure, if the 
handbook data are used. We estimate it to be less than 
0.3 % for separations larger than 50 nm. 



VI. APPLICATION TO OPTICAL DATA OF 
REF. (2^ 

We apply now the windowed dispersion relation Eq. 
([Mil to some of the optical data for gold films of Ref. 
[23 |. As we discussed at the beginning of Sec. IV these 
data are particularly interesting for our purposes, be- 
cause the window method is expected to provide very re- 
liable results for the dielectric permittivity e(«^) of these 
samples, thanks to the fact that in [22| both the real and 
the imaginary parts of e(a;) were measured independently 
by ellipsometry over the wide range of frequencies from 



TABLE IL Values of the Drude parameters for samples 1, 4 
and 5 of Ref. 22]. 



Sample 


ujp [eV /%] 


7 [meV/?i] 


1 


6.82 ± 0.08 


40.5 ± 2.1 


4 


8.00 ± 0.16 


35.7 ± 5.1 


5 


8.38 ± 0.08 


37.1 ± 1.9 



35 meV/h to about 9 eV/h, which includes longer wave- 
lengths than the tabulated data. Of the five gold sam- 
ples described in [2^ . we considered samples 1, 4 and 5. 
Samples 1 and 4 consisted of gold films deposited on a Si 
substrate, with thicknesses of 400 nm and 120 nm respec- 
tively, while sample 5 consisted of an annealed gold film 
with a thickness of 120 nm deposited on a cleaved mica 
substrate. The Drude parameters for these films were 
determined in ^23\ by several different methods. For the 
convenience of the reader they are reported in Table |lT] 
We note that samples 1 and 5 (see Table II of [l^]) are 
characterized, respectively, by having the smallest and 
largest values of the plasma frequency iOp, and therefore 
the corresponding dielectric functions display the 

largest and the least deviation from the handbook data 
(see Fig. 11 of Ref. 



2211). Sample 4 on the other hand is 
the sample for which the Drude parameters were found 
to have the maximum uncertainties. 

We evaluated the dielectric functions e{i(,) for these 
three samples, by using our windowed dispersion formula, 
and the traditional KK formula. In Fig. [T3] we display 
the fractional difference (ewin(«C) ~ eKK(j?))/eKK(«C) (in 
per cent) between the windowed and the KK estimates 
of the dielectric functions for these three samples, versus 
the imaginary frequency. When computing, by either 
method, the quantity eexpt(*C) in Eq. ([T8| . the hand- 
book data were used for frequencies larger than 9 eY/h. 
This is not expected to have much of an impact on the 
obtained results, because these large frequencies do not 
contribute much to e(i^) in the considered range of imag- 
inary frequencies. We observe first of all that the win- 
dow estimates for e(z^) corresponding to the two choices 
of the parameter b agree within half per cent with each 
other. This fact proves the good degree of KK consis- 
tency of the optical data of [22|. Interestingly, however, 
the three samples display significant differences for what 
concerns agreement of the window and KK estimates. 
While for sample 5 the window estimates agree within 
less than one per cent with the KK values, for sample 1 
and sample 4 significant deviations up to 6 % and 10 % 
respectively are found, the window approach providing 
systematically larger values than the KK method. By 
looking at the thick solid lines in Fig. [SJ one can see that 
the obtained discrepancies are too large to be explained 
by uncertainties in the Drude parameters for these film, 
listed in Table [IlJ We speculate that the found discrepan- 
cies between the windowed and the KK values obtained 
for samples 1 and 4 are related with deviations of the 
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FIG. 14: Per cent difFerence between the KK and the window 
estimates for the permittivity e(i^) for samples 1, 4 and 5 of 
Ref. [121. The window function is as in Eq. p3p . Solid lines 
are for 6=1 eV/h and dashed lines for b — 1.5 eV/h. The 
Drude parameters used to extrapolate the data are listed in 
Table [ni 

dielectric function for these films from the Drude model, 
evidenced by the presence reported in [22|] of unexplained 
absorption bands in the IR, for the samples (1,2,3 and 
4) deposited on Si substrates (see Fig. 7 of [22]). The 
large errors in the Drude parameters that were obtained 
for sample 4 probably reflect the particular inadequacy 
of the Drude model for this sample. It is quite possible 
therefore that the KK estimate of e(i^), which heavily re- 
lies on the Drude extrapolation (see Fig. 01), is much less 
accurate for samples 1 and 4. In Fig. [15] we display the 
fractional difference [Pwin(a,T) — PKK{a,T)]/ PKK{a,T) 
(in per cent) between the windowed and the KK pres- 
sures versus plate separation (in microns) for samples 1, 
4 and 5. As we see the window and the standard KK pre- 
dictions for the Casimir pressure agree within less than 
0.1 per cent in the case of sample five, while for sample 1 
and for sample 4 the window method predicts a stronger 
attraction than the standard KK approach. The theo- 
retical uncertainty in the window prediction of the pres- 
sure, that can be estimated by considering the difference 
between the obtained pressures for the two considered 
values of the window parameter &, is less than 0.15 per 
cent for all samples and all displayed separations. 



VII. PLASMA PRESCRIPTION 

In recent years it has been claimed that the Drude 
prescription leads to results that are in contradiction 
with precise measurements of the Casimir force in the 
submicron region [ll|. However, the recent experiment 
[13} has been interpreted as being fully consistent with 
the Drude prescription. Independently of its accordance 
with experiments it has been claimed that the Drude 
prescription, when used at cryogenic temperatures in the 
idealized case of samples without defects, leads to a vi- 
olation of the Nernst heat theorem |8l|. For a detailed 
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FIG. 15: Per cent difference between the window and the KK 
estimates for the Casimir pressure between two gold plates, 
for samples 1, 4 and 5 of Ref. The windowed pressure 

Pwin was computed using Eq. (fT3)) . with 6 = 1 eV/h (solid 
line) and b = 1.5 eV/h (dashed line). The Drude parameters 
used to extrapolate the data are listed in Table [III 

and updated discussion of these topics we address the 
reader to the monograph 2]. It has been argued that 
the above difficulties point to a fundamental flaw in the 
Drude approach, connected with relaxation processes of 
conduction electrons, and it has been suggested that the 
above deficiencies of the Drude prescription can be re- 
solved if relaxation processes of conduction electrons are 
neglected when evaluating Lifshitz formula. According 
to this new scheme, the quantity e(z^) to be plugged in 
Lifshitz formula should not be identified with the actual 
dielectric function of the metal, which includes such dissi- 
pative processes. One should rather use a modification of 
the dielectric function, so constructed as to disregard re- 
laxation processes of conduction electrons. The modified 
expression of the permittivity, given in Eqs. and (|25)) 
below, has been dubbed generalized plasma model, and 
following [2,] we shall denote it by egp(if ). The name orig- 
inates from the fact that within this model conduction 
electrons are described by the undamped plasma model 
ep(a;) of infrared physics 

ep{u;) = 1 - ^ , (22) 

to which the Drude model Eq. ([2]) reduces when the dis- 
sipation parameter 7 is set to zero. We note that the 
generalized plasma model is not completely immune of 
troubles either, because its 1/cj^ singularity entails viola- 
tion of the Bohr- van Leeuwen theorem of classical statis- 
tical physics, which is expected to hold for normal metals 
[To| . In Ref. [2| the following formula was obtained for 
egp(«0: 

egp(»0 = 1 + tI + - / ^-7^ 4orc(^) , (23) 

? Jo + <; 

where the second term on the r.h.s. accounts for conduc- 
tion electrons, now described as a collisionless plasma. 
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and e'^^j.^{uj) denotes the contribution of core electrons. 
The plasma frequency occurring in Eq. must be ex- 
tracted from the data. In the case of gold, the authors 
of TT] adopted the values ujp = 8.9 eV/U and 7 — 0.0357 
eV/h, that were obtained from resistivity measurements 
performed at several different temperatures on their sam- 
ples. 

An equivalent representation of the generalized plasma 
dielectric function egp(«^), using the windowed dispersion 
relation Eq. p^ . can be obtained by using the following 
representation for ecorc(^^): 



M) 




duj- 



duj- 



62 



(24) 



The equivalent representation of egp{i^) then follows if 
we replace the second term on the r.h.s. of Eq. (^5)) 
representing the core contribution, by the above repre- 
sentation of ecore(«C)- 



— / du!- 



duj 



(25) 



If the frequency b is chosen to be smaller than the on- 
set frequency Winter of interband transition, the quantity 
e^Qj.(,(a;) can be regraded as a constant Ccorc- The first 
integral on the r.h.s. of the above Equation can then be 
evaluated exactly. After doing it we obtain the following 
formula for egp(if): 





dij- 



(26) 



The above formula can be used as a substitute for Eq. 
(|23p to compute egp(i^) starting from optical data. 



VIII. CONCLUSIONS AND DISCUSSION 

In recent years, much attention has been devoted 
to theoretical and experimental studies on the thermal 
Casimir effect between two metallic bodies. From the 
theory point of view the problem is that of understanding 
whether when computing the Casimir force, conduction 



electrons should be modelled as a collisionless plasma, as 
it is done in infrared physics, or whether ohmic relaxation 
processes should be fully considered. The two alterna- 
tives have been dubbed as plasma and Drude prescrip- 
tions, respectively. Depending on the chosen approach, 
different predictions result for the magnitude of the ther- 
mal Casimir force. The present experimental status of 
the problem is perplexing, because contradictory results 
have been obtained by different experiments. On one 
hand, a series of experiments (Tl| using a micro-torsional 
oscillator to measure the Casimir force in the short sep- 
aration range from 160 nm to 750 nm, were shown to be 
in agreement with the plasma prescription and to rule 
out the Drude prescription. Also a large distance torsion 
balance experiment [12] in which the Casimir force was 
measured in the separation range from 0.48 to 6.5 /im, 
obtained results that are in agreement with the assump- 
tion of ideal metal plates, and are in contradiction with 
the Drude model. On the contrary a new large distance 
torsional balance experiment [l3| . probing the Casimir 
force between a large spherical lens and a flat plate in 
the range from 700 nm to 7 fim, was claimed to be fully 
consistent with the Drude prescription, and to rule out 
the plasma model. 

Motivated by the puzzling problem of the thermal 
Casimir effect, in the present paper we have developed 
mathematical tools that permit to obtain reliable predic- 
tions for the Casimir force in experiments using metal- 
lic plates, like those described above. This is an indis- 
pensable requisite for an assessment of the plasma versus 
Drude conundrum, when sub micron separations are con- 
sidered. This is so because for separations smaller than 
half a micron the plasma and the Drude models predict 
magnitudes of the Casimir force that differ by just a few 
per cent, and therefore it is essential to make sure that 
the Casimir force can be actually predicted with per cent 
precision or better. Our work aims at resolving a diffi- 
culty that arises when one tries to obtain an accurate pre- 
diction for the Casimir force, starting from optical data 
for the plates of the Casimir apparatus. We considered 
specifically the case of gold films, which are used in the 
experiments described above. As it is well known, in or- 
der to compute the Casimir force, one needs to know the 
dielectric function e(i^) of the involved materials, along 
the imaginary frequency axis. This function cannot be 
measured directly, and usually it is computed by means 
of a KK formula, which expresses e(i^) in terms of optical 
data for the (imaginary part of the) measurable dielec- 
tric function e{uj). The difficulty addressed in this paper 
originates from the fact that optical data of gold are as a 
rule available only for wavelengths A shorter than a max- 
imum wavelength belonging to the IR region. This is 
the case for example for the widely used handbook data 
in [2Qj, where the complex index of refraction for gold 
is tabulated for wavelengths smaller than about 10 mi- 
crons. The same is also true for the recent measurements 
of [22I . where the dielectric function e(aj) for several gold 
films was measured by ellipsomtery in the range from 33 
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to 0.14 fj,m. The problem is that knowing the dielectric 
function of gold for these wavelengths, one cannot obtain 
an accurate estimate of the KK integral for e(i^). This 
so because the 1/uj singularity displayed by the dielec- 
tric function of ohmic conductors, makes low frequencies 
give a very large contribution to the the KK integral. Un- 
til now this difficulty has been resolved by extrapolating 
available optical data towards zero frequency by means 
of the simple Drude formula, which is known to provide 
a reasonable approximation to the dielectric function of 
ohmic conductors below the interband transition. Since 
the Drude extrapolation gives a very large contribution to 
the KK integral, the obtained predictions for the Casimir 
force are affected by a significant uncertainty, caused by 
experimental errors in the values of the Drude parame- 
ters. It was estimated in [2^ that experimental uncer- 
tainties in the Drude parameters easily determine a one 
per cent uncertainty in the Casimir force between two 
gold plates at separations about 100 nm. One cannot 
hope to resolve the problem by further extending opti- 
cal measurements to lower frequencies, because for this 
purpose one would have to reach the Tera-Hertz region, 
which is practically very difficult. 

Elaborating on an earlier proposal by the author psf . 
we have shown in this paper that it is possible to resolve 
this difficulty by using modified dispersion relations to 
evaluate e(iC). The new form of dispersion relations in- 
volves certain weight functions, that we called window 
functions, that are aimed at suppressing the contribu- 
tion of low frequencies, for which no optical data are 
available. The key feature of the new dispersion rela- 
tions, as contrasted with the standard KK formula, is 
that they involve both the real and the imaginary parts 
of the permittivity. This does not constitute a problem, 
because refined optical techniques like ellipsometry per- 
mit to measure both quantities accurately [22| . The win- 
dow functions introduced in our first work produced un- 
satisfactory results, when applied to the handbook data 
for gold [15j. In this paper, we demonstrated by means 
of a Monte Carlo simulation that the bad performance 
of these window functions originated from the fact that 
they lead to a large amplification of possible small exper- 
imental errors in the optical data. In this paper we have 
addressed this instability by introducing a new class of 
weight functions, which differently from the old ones lead 
to a positive definite kernel in the dispersion formula, in 
such a way that no instability occurs anymore. We have 
shown that the improved choice of window functions effi- 
ciently suppresses the contribution of low frequencies, in 
such a way that it is now possible to compute with high 
precision the Casimir pressure between two gold plates in 
vacuum, using currently available optical data, like the 
handbook data [13] and the new data of Ref. . 

When we used the handbook data, we found it neces- 
sary to exclude from our computations all data for fre- 
quencies less than 0.6 eY/h, whose original source was 
Ref. jl^ . The reason for not using these IR data is that 
the film used by these authors appears to posses a plasma 



frequency of 7.7 eV/Ti, which is exceedingly smaller than 
the value of 9.2 eV/h that we obtained from Theye's 
data [lOl , which are part of the handbook data from the 
near IR to the UV. The latter value of the plasma fre- 
quency is very close to the commonly accepted value of 
9 eY/h [24]. After excluding the data of [23] the window 
method produces less accurate results, because the resid- 
ual handbook data that ar left over start from a rather 
large frequency of 0.6 eV/h. Nevertheless, the Casimir 
pressure can still be estimated by the window method 
with an uncertainty that we estimated smaller than 0.3 
%. When the Drude prescription is used, the values of the 
Casimir pressure obtained by the window method agree 
with those obtained by the standard approach, based on 
KK relations and the Drude extrapolation, within 0.2 
% for separations larger than 50 nm, and therefore they 
remain inconsistent with the experimental data of llj. 

A much higher precision can be obtained when the data 
of [22j] are used. Thanks to the fact that these data ex- 
tend to longer wavelengths than the tabulated data, the 
uncertainty in the windowed prediction for the Casimir 
pressure decreases to 0.1 %, for separations larger than 
50 nm. This is an order of magnitude better than the 1 % 
uncertainty that was estimated in [22| to result from un- 
certainty in the Drude parameters, within the standard 
KK approach. The magnitudes of the Casimir pressure 
predicted by the window method were found to be in very 
good agreement with predictions based on the standard 
KK formula for the mica sample of . However for the 
two samples deposited on Si substrates that we consid- 
ered significant discrepancies were found, the windowed 
Casimir pressures exceeding the KK values by more than 
1 % for sample 1 and by more than 2 % for sample 4, for 
separations less than 200 nm. The discrepancy for sam- 
ples 1 and 4 is too large to be explained by uncertainty 
in the Drude parameters, and we attribute it to the fact 
that the Drude model was shown in 22| not to describe 
well the dielectric function of these films in the IR, as 
evidenced by the presence in their optical data of an un- 
explained IR absorption band. It is quite possible that 
for this reason predictions for the Casimir force based on 
the Drude model are less accurate for these films. 

In conclusion our results demonstrate that the window 
approach can be used to predict the Casimir pressure 
between two metallic plates very accurately, on the basis 
of standard optical data extending from the IR to the 
UV, without the need of further extending optical mea- 
surements to longer wavelengths than those reached by 
current laboratory optical apparatus. We hope that the 
results of this work will help clarifying present difficulties 
raised by the interpretation of experiments probing the 
thermal Casimir effect. 
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